Possible High-Energy Neutrinos from the Cosmic Accelerator 

RX J1713. 7-3946 



J. Alvarez-Muniz 

Bartol Research Institute, 
University of Delaware, Newark, DE 19716, USA. 

and 

F. Halzen 

University of Wisconsin, Department of Physics, 
1150 University Avenue, Madison, WI 53706, USA. 

ABSTRACT 

The observation of TeV-gamma rays of 7r° origin from the supernova remnant 
RX J1713. 7-3946 might have revealed the first specific site where protons are 
accelerated to energies typical of the main component of the cosmic rays. In 
this letter we calculate the high-energy neutrino flux associated with this source 
to be 40 muon-type neutrinos per kilometer-squared per year. We perform the 
same calculations for other known sources of TeV-gamma rays and show how 
neutrino observations can establish whether the TeV-gamma rays emitted by 
blazars and supernova remnants are the decay products of neutral pions and 
thus unequivocally establish the sources as cosmic accelerators. 

Subject headings: acceleration of particles — cosmic rays — gamma rays: theory 
— neutrinos — supernova remnants 

1. Introduction 

Supernova remnants (SNRs), such as the Crab (Weekes et al. 1989), Cassiopeia A 
(Aharonian et al. 2001) and blazars, such as Mkn 421 (Punch et al. 1992) and Mkn 501 
(Quinn et al. 1996), emit gamma rays with energies up to tens of TeV, (see (Ong 1998), 
(Catanese & Weekes 1999), and (Pohl 2001) for reviews). Unfortunately, these observations 
cannot by themselves establish whether these high-energy photons are produced by electron 
beams via inverse Compton scattering (Pohl 2001), or by proton beams via the production 
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and subsequent decay of neutral pions (Gaisser et al. 1995; Pohl 2001). The mechanisms 
for the production of the highest energy photons by purely electromagnetic processes do not 
naturally accommodate TcV-photons but can be stretched to explain the observed spectrum. 
The question is central to the problem of the origin of cosmic rays and it has been pointed 
out that neutrino observations can unequivocally settle this debate (Gaisser et al. 1995). 
In hadronic mechanisms the high-energy photons are accompanied by neutrinos because 
charged pions are produced along with 7r°'s. 

Recently the CANGAROO collaboration has produced evidence for the observation of 7- 
rays of tt^ origin from the SNR RX J1713. 7-3946 (Enomoto et al. 2002). We here show how a 
lower limit on the neutrino flux from this source can be established from energy conservation 
and routine particle physics. We show that this limit is accessible by a northern hemisphere 
neutrino observatory such as Antares (Feinstein 1999) operated over several years. Obvi- 
ously, if photons are partially absorbed in the source or on the interstellar medium, this 
calculation underestimates the actual neutrino flux. 

The true nature of RX J1713. 7-3946 is still a matter of discussion (Butt et al. 2001). 
Our calculation illustrates the possibility to disentangle electromagnetic and hadronic sources 
of high-energy radiation observationally. We will establish, in a largely model-independent 
way, a direct relation between the lower limit on the neutrino flux and the TeV-gamma 

ray observations for sources other than RX J1713. 7-3946. We compute the response of 
a kilometer-scale neutrino observatory such as IceCube (IceCube 2002) to that flux, and 
conclude that the event rates are adequate to confirm or rule out the hadronic origin of 
the highest energy photons. Although the flux from a single source may still be small, this 
conclusion is credible because a neutrino telescope will be operated for a decade. 

We should emphasize that, even if RX J1713. 7-3946 is a truly cosmic ray accelerator, 
it is likely that the dominant source of the cosmic rays has not been identified. The higher 
energies might be associated with older supernova remnants that may be astronomically 
uninteresting and have escaped attention. Because of their restricted angular resolution 
atmospheric Cherenkov telescopes have scrutinized a limited number of sources. Moreover, 
TeV-gamma ray instruments are limited to nearby sources because of absorption of the 
flux on interstellar light. A Mediterranean telescope (Feinstein 1999; Gricder et al. 2001; 
PiateUi 2002) and AMANDA/IceCube (Andres et al. 2000; IceCube 2002) wiU survey the 
whole sky. 
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2. Neutrino flux from gamma ray observations 

Let's set out to test the assumption that TeV photons produced by a supernova remnant 
or by an active galaxy are of hadronic origin, i.e. that they are the decay products of 
7r°'s. Neutral pions can be produced by accelerated protons in either pp or pry interactions 
depending on the relative target density of photons and protons in the source region where 
the protons are accelerated (Halzen & Hooper 2002). The + neutrino flux [dNj^/dEi) 
produced by the decay of charged pions in the source can be derived from the observed 7-ray 
flux (dN-y/dEj) by imposing energy conservation (Stecker 1979; Stecker & Salamon 1996): 



/■^r'' dN Z-^^" dN 

/ E,^dE,^K E^^dK, (1) 

where E"^^"^ {E^^^) is the minimum (maximum) energy of the photons that have a hadronic 
origin. E™™ and i^™^^ are the corresponding minimum and maximum energies of the neu- 
trinos. The factor K depends on whether the tt^'s are of pp or p^ origin. Its value can be 
obtained from routine particle physics. In pp interactions 1/3 of the proton energy goes into 
each pion flavor on average. In the pion-to-muon-to-electron decay chain 2 muon-neutrinos 
are produced with energy Et^/A for every photon with energy Et,/2 (on average). Therefore 
the energy in neutrinos matches the energy in photons and K — 1. 

Important parameters in relating the neutrino flux to the 7-ray flux are the maximum 
and minimum energies of the produced photons and neutrinos. They appear as limits in 
the integrals in Eq. (1). The maximum neutrino energy is fixed by the maximum energy 
of the accelerated protons (E^^) which can be conservatively obtained from the maximum 
observed 7-ray energy E^^. Following our previous assumptions: 

1 2 

for the pp case. The minimum neutrino energy is fixed by the threshold for pion production. 
For the pp case: 

(2mjj + TTT-Tr)^ — 2ml 
=^ ' 

where F is the Lorentz factor of the accelerator relative to the observer. One has to keep 
in mind that in blazar jets the shocked regions where protons can be accelerated move 
relativistically along the jet and are therefore boosted relative to the observer. The average 
minimum neutrino energy is obtained from using Eq. (2). 

For p7 interactions it can easily be shown, using the A-resonance approximation, that 
K — A. The minimum proton energy is given by the threshold for production of pions 
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via the A-resonance, and the maximum neutrino energy is 5% of the maximum energy 
to which protons are accelerated. More detailed calculations will be published elsewhere 
( Alvarez- Muniz & Halzen 2002). 



2.1. Neutrino events from RX J1713.7-3946 

Recently the CANGAROO collaboration has shown that the energy spectrum of the 
7-ray emission from the supernova remnant RX J1713. 7-3946 matches that expected if the 
7-rays are the products of vr^ decay generated in pp collisions. We normalize the neutrino flux 
from this source using Eq. (1). We calculate the total energy emitted in 7-rays integrating 
the 7-ray spectrum from vr" decay calculated by (Enomoto ct al. 2002). Gamma rays of 
energies up to ~ 10 TcV have been observed, indicating that protons are accelerated up to 
^max ^ rj^^Y ^2) then gives E^"^^ ~ 10 TeV. Ef'' is obtained from Eqs. (2) and (3), 
where we have taken F ~ 1 assuming the velocity of expansion of the supernova shell is not 
much larger than ~ 10,000 km/s. The data in (Enomoto et al. 2002) has been fitted with a 
model that uses a ~ E~'^ input proton spectrum. We assume that all charged pions decay in 
the environment of the supernova and hence the neutrino spectrum follows the input proton 
spectrum. Under these assumptions our calculated + 9^ flux is given by: 



= 4.14 X 10-" (E^/l TeV)-' [cm^^ s"' TeV"']. (4) 



The neutrino event rate is calculated in a straightforward manner by convoluting the 
Ufj_ + z/^ flux with the probability of detecting a muon produced in a muon-neutrino interaction 
(Gaisser et al. 1995): 

A^events - Aff X T X J '^{E,)P,^^{E,)dE,, (5) 

where T is the observation time and A^f^ is the effective area of the detector. 

The neutrino event rate in a detector of effective area = 1 km^ is shown in Fig. 1 
as a function of the muon energy threshold of the detector [E^^^). Also shown is the at- 
mospheric neutrino background in a 1° x 1° window in the sky. The dashed and dotted 
curves show the contribution from vertical downgoing and horizontal atmospheric neutri- 
nos respectively, bracketing the maximum and minimum values of the expected atmospheric 
neutrino background (Volkova 1980). The actual background depends on the (in general) 
variable zenith angle of the source as seen from the detector. Clearly the background does 
not pose a problem for the identification of this source. Due to the much larger muon flux 
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produced in cosmic ray interactions in the atmosphere compared to the expected muon flux 
from neutrino interactions, only neutrinos going through the Earth can be identified. For 
this reason a southern hemisphere source such as RX J1713. 7-3946 is only accesible to a 
northern hemisphere neutrino telescope. 

In a few years of operation a detector in the northern hemisphere such as the current 

Antares {A^e ~ 0.1 km^, E^^' ~ 40 GcV) will be able to identify RX J1713.7-3946 with a 
Poisson chance probability that the event rate is a fluctuation of the atmospheric neutrino 
background smaller than 1%. The IceCubc neutrino telescope {E^^^^ ~ 200 GeV) will do 
the same job in just 1 year of operation if a source similar to RX J1713. 7-3946 exists in the 
northern hemisphere. Besides, these detectors have the potential to discover new SNR's that 
might have escaped the scrutinity of 7-ray Cherenkov telescopes. 

It should be mentioned that we have not attempted to correct for absorption of the 7-ray 
flux in the source or in the interstellar medium. The latter is justified since the distance to 
the source is poorly known. For these reasons our calculation could well be conservative. We 
have not included the possible effects of oscillation of the cosmic beam (Crocker et al. 2000). 

The claims made by the CANGAROO collaboration have been challenged in a recent 
paper (Reimer et al. 2002). The statement is that the EGRET experiment detected less 
gamma rays from the associated source 3EG J1714-3857 than expected from 7r° decay. How- 
ever, the EGRET measurements were made between 1991 and 1995 (Hartman et al. 1999), 
i.e. several years before the CANGAROO observations which were carried out in 2000-2001. 
It is possible that the source might have experimented large flux variations over the period 
of time separating both observations. It is not the purpose of om work to further discuss this 
point. We stress here that neutrino telescopes have the sensitivity to settle this important 
debate. 



2.2. Are there other proton cosmic accelerators? 

Gamma rays at TeV energies have also been detected from other shell-type SNRs such 
as the Crab Nebula, Cassiopeia A and upper limits have been reported from the SNR-like 
shell Sagittarius A East at the Galactic center (Buckley et al. 1997). Besides, TcV 7-rays 
have been observed from at least two blazars: Mkn 421 and Mkn 501. For all these objects 
the observations are less constrained than in the case of RX J1713. 7-3946, and in fact models 
in which the 7-rays are produced by inverse Compton scattering of ultra-high energy elec- 
trons on ambient photons can be stretched to explain the TeV observations. For Cassiopeia 
A and Sagittarius A East we proceed as in the case of RX J1713. 7-3946 normahzing the 
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neutrino flux by the 7-ray flux expected from 7r° decay calculated in (Ellison et al. 1999) 
and (Melia et al. 1998) respectively. For the Crab we assume that the measured high energy 
7-ray spectrum has been produced in pp collisions, and the corresponding gamma flux from 
the Mkn 421 and Mkn 501 in p'-f collisions, and we compute the accompanying neutrino flux 
and neutrino event rate using Eq. (1). 

We perform the actual calculations by assuming that both the photon and the neutrino 
fluxes are generated by protons with a differential energy spectrum proportional to E^'^, 
and assume the neutrino spectrum is parallel to the proton spectrum. The spectral index of 
photons is taken from the observations assuming the 7-rays are of hadronic origin. 

In the latter calculations we, conservatively, only identify the observed 100 GeV-TeV 
7-rays with 7r° origin. This spectrum does of course continue down to low energies and turns 
over at m^o/2, as it does for RX J1713. 7-3946. This low energy spectrum is not revealed 
by the observations and modeling it is not straightforward because of possible cascading, 
boost factors, etc. We therefore prefer to limit the spectral information to the region of 
observation, (corresponding to photon energies above a certain value which we denote by 
E^^^), thus producing a lower limit on the neutrino flux. 

High-energy photons can be absorbed in the interstellar infrared radiation. Because 
the amount of absorption depends on the distance to the source and on the poorly known 
infrared background (Stecker & de Jager 1998) we are simply going to use as the max- 
imum observed energy keeping in mind that this assumption leads to an underestimate of 
the neutrino flux. We have used the measured TeV 7-ray fluxes given in the following 
references: Crab (Aharonian et al. 2000), Mkn 421 (Krennrich et al. 2001), and Mkn 501 
(Aharonian et al. 1999). 

The event rates per km^, and the values of the relevant parameters are shown in Table 1. 
We calculated the event rate for two choices of E^^*, namely, the energy threshold of the 
Cherenkov telescope that performed the measurement, and E^'^'^ = 100 GeV, illustrating the 
sensitivity of the neutrino event rate to E^'^^. Also tabulated are the background of atmo- 
spheric neutrinos above i^™™ in a 1° x 1° pixel in the direction of the corresponding source, 
and the Poisson chance probability of obtaining a background of events from atmospheric 
neutrinos equal to the total expected signal during the time of the observation. 

For the SNR Cassiopeia A we obtain 5 events per km^ per year with a 4.7% chance 
probability that the rate is a fluctuation of the atmospheric background. Due to the extreme 
inferred energetics of Sagittarius A East, about 40 times the energy released in a typical 
supernova explosion (Meha et al. 1998), the neutrino event rate is large: ~ 140 km~^ yr~^. 
Sagittarius A East is however too close to the Galactic center to be observed by a detector 
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at the South Pole but it could be detected in a northern hemisphere telescope. 

We conclude that the hadronic origin of the observed TeV-gamma rays can be con- 
firmed or disproven with confidence given that neutrino telescopes will be operated over 
many years. Because of the total transparency of the Universe to neutrinos, we also expect 
a larger number of sources. RX J1713. 7-3946 awaits to be detected by northern hemisphere 
telescopes (Feinstein 1999; PiateUi 2002; Balkanov 2000; Grieder et al. 2001). Its detec- 
tion will be a milestone in neutrino astronomy. Our results are also encouraging for the 
AMANDA/lceCubc detectors at the South Pole. They clearly have the capability to iden- 
tify cosmic accelerators observed in 7-rays. They may also reveal sources that have escaped 
detection by Cherenkov telescopes. 

We thank Diego F. Torres for discussions. We are grateful to the anonymous referee 
for valuable comments and suggestions on the manuscript. J.A.-M. is supported by NASA 
Grant NAG5-10919 and thanks the Department of Physics, University of Wisconsin, Madison 
where part of this work was done. This work was supported in part by DOE grant No. DE- 
FG02-95ER40896 and in part by the Wisconsin Alumni Research Foundation. 
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Fig. 1. — Neutrino event rate per year from RX J1713. 7-3946 as a function of the muon 
energy threshold in a km-scale neutrino telescope. Also shown is the event rate expected in 
a 1° X 1° bin in the sky from atmospheric neutrinos. The dashed curve shows the predic- 
tion for vertical downgoing atmospheric neutrinos and the dotted the corresponding one for 
horizontal atmospheric neutrinos. 
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Tabic 1: N^^^^y^ events per km^ from the Crab and Mkn421 in one year and from Mkn 501 
during the 8 months of flaring activity in 1997. We are assuming that protons of energy up 
to 10^° eV can be accelerated in the Mkn's. Also shown is the expected rate of atmospheric 
neutrino events (assuming Volkova's flux) in a 1° x 1° bin on the sky around the position 
of the source (Macomb & Gehrels 1999) as seen from IceCube. The last row is the chance 
probability of observing a background of atmospheric neutrino events equal to the total 
signal in the same period of time. See text for references. 



Source 


Crab 


Mkn 421 


Mkn 501 


[GeV] 

^ma^ [TeV] 

TV ^- 

yyatm 

Poisson prob. (%) 


100 500 
20 20 
10 10 
9.5 4.0 
3.5 3.5 
1.2 18.8 


100 260 
17 17 
10^ 10^ 
7.5 3.5 
1.5 1.5 
0.1 3.8 


100 500 
24 24 
10^ 10^ 
8.0 5.0 
1.0 1.0 
0.001 0.3 



